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Abstract. Analytic formulae are presented to construct 
detailed secular lightcurves of both early asymptotic gi- 
ant branch (AGB) and thermally pulsing AGB stars. 
They are based on an extensive grid of evolutionary cal- 
culations, performed with an updated stellar evolution 
code. Basic input parameters are the initial mass Mi, 
0.8 < Mi/M & < 7, metallicity Z x = 0.0001,0.008,0.02, 
and the mixing length theory (MLT) parameter. The for- 
mulae allow for two important effects, namely that the 
first pulses do not reach the full amplitude, and hot bot- 
tom burning (HBB) in massive stars, which are both not 
accounted for by core mass - luminosity relations of the 
usual type. 

Furthermore, the dependence of the effective temper- 
ature and a few other quantities characterizing the con- 
ditions at the base of the convective envelope, which are 
relevant for HBB, are investigated as functions of luminos- 
ity, total and core mass for different formulations of the 
convection theory applied, MLT or Canuto & Mazzitelli's 
(|1991D theory. 



Key words: Stars: AGB and post-AGB - Stars: evolu- 
tion - Stars: fundamental parameters - Methods: numer- 
ical 



1. Introduction 

As a matter of fact, full AGB stellar evolution calculations 
are unable to provide the statistical information needed for 
purposes like population synthesis. There are three main 
reasons for this. 

First, the calculations are rather lengthy, reflecting the 
complexity of the inner structure of AGB stars and of 
the temporal evolution due to thermal pulses (TPs , also 
called helium shell flashes; see e. g. Iben & Renzini 1983| ; 
Boothroyd & Sackmann 1988] , in the following BS88; Lat- 
tanzio 



1986, LA86; Wagenhuber & Weiss 1994.WW94 



and Blocker |1995| , B95) 



Second, the calculations involve at least two poorly 
known parameters, the mass loss efficiency rj and the mass 
loss prescription, and the mixing length (MLT) parameter 
a, each of them only meaningful in the framework of the 
a priori chosen description of the mass loss and convec- 
tion theory, respectively. Since it has become clear that 
AGB evolution is dominated by mass loss (Schonberner 
1979 ), numerous formalisms giving the mass loss rate as 
a function of the stellar parameters have been applied 
(e.g. Reim ers 1 975, Iben & Renzini 1983 , Mazzitelli & 
D'Ant ona b.986j Wood & Faulkner |l986j , Vassiliadis & 
Wood |1993 VW93, B95). The uncertainties were cast into 



the "efficiency" parameter, which depends on the mass 
loss prescription used, and which is actually a parame- 
ter that can be calibrated from AGB population synthesis 
models (e.g. GdJ, Groenewegen & de Jong 1994). 

In addition, the MLT parameter is related to a series 
of problems concerning (a) the correct effective tempera- 
ture scale of late-type stars derived from observations, (b) 
the opacity at low temperatures (Alexander & Ferguson 
1994), (c) the time-dependent problem of dust formation 
(Fleischer et al. 



1992 



and references therein), and (d) the 
theory of convection itself. 

The third restriction of the use of full evolutionary 
calculations is the occurrence of the third dredge-up as a 
consequence of TPs, which is considered to be responsi- 
ble for the formation of carbon stars. However, canonical 
stellar evolution calculations do not consistently predict 
the formation of carbon stars in the mass range that they 
are observed, roughly between 1.2-1.5 AIq depending on 
metallicity and about 5 M Q (e.g. GdJ, Groenewegen et 
al. 1995). The explanation is still under debate and r ecent 
progress has been made (see e. g. Frost & Lattanzio 1996 , 
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Staniero et al. 1997, Herwig et al. 1997) but not to the 
degree that the results can easily be included in synthetic 
AGB models. Therefore these models introduce a fudge 
parameter to describe the third dredge-up, in particular 
the dredge-up efficiency A. 

To summarize: In order to describe the evolution of 
whole stellar populations using synthetic calculations, one 
has to explore a parameter space with at least five dimen- 
sions: (Mj, Zi, a, rj, A). This has been done in the past e. g. 
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by Renzini fc Vo li ( |198l[ ), J0r gensen ( |1991[), Gr oenewegen 
& de Jong ( p.993| , GdJ) and Marigo et al. ( |1996[ ). Although 
these works have contributed a lot to an improved under- 
standing of AGB evolution, there are some obvious short- 
comings: e.g. the data collected by GdJ and frequently 
used since then are partially incomplete, making some ad 
hoc inter- and extrapolations necessary, since results ob- 
tained in three decades are combined with each other, and 
they arc oversimplified in some respects, e. g. concerning 
hot bottom burning (HBB). Furthermore, some data such 
as the rapid luminosity variations during a TP were ne- 
glected by many authors, although they are needed for 
some of the applications mentioned below. 

The aim of this paper is to provide essential theoretical 
data for applications like classical synthetic AGB evolu- 
tion. In this sense this paper provides a fully updated and 
improved set of relations with respect to GdJ. 

The paper is organized as follows. In the next sec- 
tion some terminology is introduced, and the intrinsic er- 
rors due to physical assumptions are discussed. After this, 
the full evolution calculations up to and along the early 
(E-)AGB are shortly described (Sect. 3). The fourth sec- 
tion contains the recipes necessary to construct the secu- 
lar lightcurves on the TP- AGB as a function of time for 
all relevant initial total masses and metallicities. Finally, 
properties of the effective temperature and the conditions 
at the bottom of the convective envelope for various model 
assumptions are discussed in Sect. 5. 

2. General remarks and computational details 

2.1. Some definitions 

Let us first define some quantities and index labels that 
are used throughout the rest of this paper. 

A thermal pulse cycle (TPC) is the time interval from 
a local maximum of helium burning, through quiescent 
hydrogen burning, up to the next TP. All quantities that 
refer to a TPC as a whole are defined to be functions of 
the conditions at the beginning of the TPC. The first TP 
is the one in which the maximum (integrated) luminosity 
produced by He-burning exceeds the maximum H-burning 
luminosity prior to this for the first time. Previously there 
may be "pre-pulses". 

Uppercase letters denote quantities that depend lin- 
early on parameters. Time intervals (in years) are denoted 
by t, in particular the duration of a TPC, the interpulse 
time, by Tj p . Luminosity and mass (L, M) are given in so- 
lar units as usual, temperatures T in K, the abundances 
of hydrogen X, helium Y and metals Z in relative mass 
units (X + Y + Z = 1), and all others in cgs units. 

Lowercase letters are used for decadic logarithms of 
the above quantities, e. g. I = lg L/Lq, except for z = 
lg Z/0.02 being the logarithmic metallicity scaled to solar. 

The core is defined to be the part of a star inside the 
location where the local hydrogen content reaches half the 



photospheric value. Since the H-burning shell is extremely 
narrow (both in radius and relative mass units), this es- 
sentially coincides with all other definitions of the core 
mass in the literature. The core growth AM C is defined as 
/ max{^,0}dr. Usually, AM C = M c (r) - M C)0 , where 
M Ci o is the core mass at the first TP. 

The mantle is the part between the photosphere and 
the core. In terms of the relative mass content, it is almost 
identical with what is frequently called "(convective) en- 
velope" . 

One has to distinguish three different TP- AGB phases, 
defined in WW94. There are: 

(i) The transition phase from the E-AGB to the first TP, 
in which pre-pulses may take place. 

(ii) The turn-on phase, comprising about the first ten TPs. 
All global quantities are aiming at their asymptotic 
behaviour, but still the deviations are significant (up 
to 60%). This is called turn-on effect or TOE here. 

(iii) The asymptotic phase, where the global quantities have 
reached their asymptotic behaviour. 

The MLT parameter a may be defined either for the 
MLT or Canuto & Mazzitelli's ( |l99f| ) theory (CMT) and 
is labelled accordingly, if necessary. 

The three sets of calculations for Z\ = 0.02, 0.008 and 
10~ 4 are called pop I, LMC and pop II in the following, 
respectively. 

Subscripts: 

b: Quantities defined at the bottom of the convective mantle, 
c: concerning the core and 
m: the mantle. 

A,B,C,D: Denote luminosity extrema in the TPC (Fig. [I]): 
A - the "slow maximum" during quiescent H-burning, 
B - the "rapid dip" following a TP, 
C - the "rapid peak" after this, 

D - the "slow dip" at the transition from He- to H-burning. 

i: Initial values on the ZAMS, and 

f: final values after the TP- AGB phase. 

H: Something produced by or related to hydrogen or 

He: helium burning. 

0: Quantities at the beginning of an evolutionary phase 

as a whole, like the E-AGB and the TP- AGB, 
TP:at the beginning of a TPC, and 
1: at the end of the corresponding phase. 
*: Denotes "effective" values measured at the photosphere. 

2.2. Intrinsic errors and physical assumptions 

About 40 model sequences were run, with varying ini- 
tial mass, metallicity or physical input parameters. Ini- 
tial masses are between 0.8 < Mi < 7, the metallicity is 
Zi = 0.02, 0.008, 0.0001. The standard values for Y on 
the main sequence were 0.28, 0.2651 and 0.25 for pop I, 
LMC and pop II, respectiv ely. T he metal mix is accord- 
ing to Anders & Grevesse (1989) for pop I, scaled by the 
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Fig. 1. A typical TPC as produced by the code described in the 



text. The labels illustrate some of the definitions in sect. 2.1 



The time axis has been stretched around the first TP shown. 

appropriate factor for pop II. For the LMC mixture, the 
abundances of C,N and O, which are considered explicitly 
in the equation of state and the nuclear reactions, deviate 



from the solar values following Russell & Bessell (1989) 
and Russell & Dopita ( 1990 ). Opacity tables were always 
used, according to their X,Y and Z values, neglecting the 
internal metal ratios. 

During the full stellar evolution model calculations 
new input physics became available which allowed us to 
test for any differences indiced by them. Two different 
physical models were used for the pop I calculations and 
Mi € {1, 1.5, 2, 3, 5, 7}: one with nuclear equilibrium rates 
and opacities from the collection of older opacity tables by 



Weiss et al. (1990), as explained in WW94, and the other 



one with a nuclear network, and more recent opacities 
taken from Rogers & Iglesias ( 1992 ), in the following re- 
ferred to as models A and B, respectively. For T < 10 4 K, 
the OPAL tables were completed by LAOL data. Both 
sets agree well around T 10 4 K. The LMC calculations 
were carried out with model A, model B was employed 
for pop II. Our nuclear network comprises the following 
species relevant for the AGB: 1 H, 3 He, 4 He, 12 C, 13 C, 
14 N, 15 N, 16 0, 17 0, 20 Ne, 24 Mg. A leakage out of the 
CNO cycle, which is unimportant for the energy gener- 
ation rate (Mowlavi et al. 1996| ), was not taken into ac- 
count. For more details, see Weiss & Truran (1990), Weiss 
et al. ( |1996D and Wagenhuber ( |1996| ). The models were 
evolved off the ZAMS through the core helium flash (for 
Mi < 2) up to the AGB without manual interventions and 
with a high numerical resolution (see WW94). 

When comparing model A and B calculations, one 
finds as the major effect that for Mj ^ 4 the initial core 
mass M c .o as a function of Mi increases slightly faster 
for model B than for model A (see Fig. 0). The great- 
est difference of O.O8M occurs for Mj = 7 (the OPAL 
data f or M c agree with those of D'Antona & Mazzitclli 
( |1996| , DAM96) to within 0.01M S ). Accordingly, the mas- 
sive stars in model B reach somewhat higher luminosities 
than their counterparts in model A with the same Mj. 



However, on the TP-AGB luminosities and time scales, 
when considered as functions of M c and AM C , agree with 
each other to within 12%, the main contribution being 
that the TOE tends to be stronger for model B. 

The convective core of the Horizontal Branch mod- 
els grow slowly into the semi-converctive regions, which, 
due to the use of the Schwarzschild-criterion, are for- 
mally stable stable in our solutions. The growth is due 
to the smoothing of (unphysical) molecular weight steps 
in the grid-control routine of the code. However, in the 
B-sequences, the growth happened in discrete "breathing 
pulses" (Castellani et al. |1985| ). The size and occurence 
depends on both physical (mixture, opacities) and numer- 
ical details. For increased resolution in both space and 
time their extent is reduced and the cores grow in a more 
continous way. 

The standard choice for the MLT parameter was a = 
1.5, which together with the adopted opacities for low 
temperatures yielded somewhat lower (i.e. the effec- 
tive temperature following our nomenclature) when com- 
pared to results by DAM96 and B95. The code used in the 
present work needs a MLT parameter approximately 0.25 
and 0.2 larger, respectively, to reproduce their results. In 
order to fit the present sun with OPAL, a — 1.578 and 
Y = 0.2803 is needed, very close to the standard values 
of 1.50 and 0.28, giving rise to a possible shift of T* by 
far within the uncertainties of the observed AGB effective 
temperature scale. Since the dependence of all results on Z 
is the dominating one, uncertainties due to the variations 
of Y for a given Z of about ±0.03 are neglected. 

Additionally, some TPCs were repeated with the CMT 
formulation of convection instead of standard MLT. This 
was done with a mixing length A proportional to the local 
pressure scale heigth H p , i. e. A = gvcmt H p . The choice 
acMT = 0.65 reproduced the effective temperatures ob- 



tained with MLT and «mlt = 1-5 very well (see Sect. 5.2). 
Beyond this, time-dependence according to the model by 



Arnett (1969) and a turbulent pressure model according 



to CMT were invoked in some test sequences. These tests 
reveal that the various convection formalisms, too, do not 
influence luminositiy and timescales by more than 2%, 
except if HBB is operating (see Sect. 4T). The time de- 
pendent convection models give rise to two interesting ef- 
fects during the "rapid luminosity peak" : T* is cooler by 
about 100K than during the quiescent phases, and Tb of 
low-mass stars, which usually never exceeds 1.5 x 10 6 K, 
reaches peak values of 5 x 10 6 K shortly after a TP, since 
the bottom of the mantle convection zone cannot retreat 
as quickly as in the case of "instantaneous" convection. 
However, L* and the TP timescales are unaffected. 

These calculations furthermore allow to estimate the 
impact of the way the initial models were generated on 
the TP-AGB. Instead of full calculations through the core 
helium flash, often the chemical profile of the last model at 
the tip of the RGB (red giant branch), or an artificial step 
profile, are used to obtain a core helium burning model by 
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Fig. 2. The lightcurves during the first TPCs for the full evo- 
lution through the core helium flash (full line) and two cases of 
artificial initial models (see text). The time axis (in Myr since 
the onset of core helium burning) of the latter two has been 
shifted by the amount indicated so that the last but one TPs 
lie on top of each other. 



means of an explicit integration method. These three pos- 
sibilities have been carried out for a pop II, M; = 1.25M Q 
model. The luminosities at the onset of the core helium 
burning phases, which last for 62, 56 and 78 Myr, are 
79.8,81.1 and 82. 3i© for the full evolution, H-profile and 
step profile cases, respectively. Only the latter experiences 
a breathing pulse. The beginning of the TP- AGB is shown 
in Fig. ||. The sequences evolved off artificial initial models 
show an irregular behaviour during the first TPCs. This 
leads to the conclusion that the composition profile of the 
H-shell provides a rather long-lasting memory, which may 
change L» and r; p at the first TP by about 10%, and will 
therefore mildly influence the evolution of low-mass stars, 
since these experience only few TPs (GdJ). 

2. 3. Outline of the method 

Here we briefly outline the method used to derive ana- 
lytical relations from the full calculations, using the core 
mass-luminosity relation as an example. 

It has been known for a long time that two im- 
portant quantities characterizing the TP- AGB, La and 
lgTjp, approximately are linear functions of the core mass 
(Paczyriski 1970, 1975). However, many different formulae 
were publis hed (e.g. Iben & Truran |l978j (IT78), Wood & 
(WZ81), LA86, BS88). 



1981 



Zarro 

To find an appropriate description of our data, first L 
and lgr were assumed to depend linearly on M c , and each 
of the three sets of calculations (pop I, II and LMC) was 
considered separately (sometimes a quadratic dependance 
was found to be more appropriate). A random subset of 
data from TPCs in the asymptotic regime without HBB, 
distributed within 0.55 < M c < 0.80, was used to obtain 
a rough estimate for the coefficients, and to determine if 
the latter are constant, or depend linearly on Z\ or lgZ;. 
One thus obtains a first guess relation g^\M c , Z{). 

Next, the residua due to the TOE are considered, still 
ignoring HBB for the moment. They can well be described 
by an exponentially decaying function of the core growth, 



i. e. the ansatz g^ k, lg \f — g^'\ is a linear function of 
AM C , with coefficients depending on M c or M Cj o and Z x or 
lg-Zi, which are determined by a non- linear fit for a sub- 
set of all the data. To continue, the least significant terms 
in the new residual function are removed. The resulting 
formula is checked against the full available data set. One 
finds that the data points previously not taken into ac- 
count are well predicted, with deviations of the order of 
less than twice the mean error with respect to the subset 
data. 

The effects of HBB are strongly non-linear functions 
of the total mass and the MLT parameter a. One cal- 
culation with a — 2 [Mi — 5M ) was carried out, and 
data from DAM96 and B95 were additionally considered. 
A large number of analytic forms with the correct quali- 
tative behaviour were tested. 

In a final step, all coefficients are improved simulta- 
neously by applying a non-linear fitting procedure to the 
complete data set. The values of the parameters are given 
with as many digits as necessary to obtain deviations of 
less than 1% with respect to results calculated with 16 
digits. The number of digits must not be mistaken as an 
indication that the parameters are "known" to this preci- 
sion, since changing one parameter usually could be com- 
pensated by changes of another one, at the expense that 
some data points were reproduced considerably worse. 

2-4- Limits of applicability 

All the formulae presented below should be considered as 
interpolations, valid for the parameter space 0.8 ^ Mi 7 
and 10~ 4 < Zi < 0.02. Below Z-, = 10~ 4 , new physical ef- 



fects occur (see Cassisi et al. 1996) that partially invalidate 
the descriptions given below. 

With the standard choice of the MLT-parameter, a — 
1.5, almost no HBB takes place, and the resulting T* are 
close to the lower limit of the what is allowed by observa- 
tions. Hence a < 1.4 is meaningless. On the other hand, 
no data with a > 2.8 are available, so that larger values 
should be avoided too. 

3. The evolution prior to the first thermal pulse 

In synthetic AGB evolution calculations, the evolution 
prior to the first TP must be taken into account to de- 
termine the total mass left at the onset of TPs. 

3.1. Up to the core helium flash 

Low-mass stars that experience the core helium flash may 
lose a significant amoun t of m ass on the first giant branch. 
J lmenez & MacDonald (|l997| ) give a core mass-luminosity 
relation (CMLR) for the RG B, w hich in terms of the vari- 
ables introduced in section |2.l| (m c — lgM c etc.) reads: 



/*=4.826+0.066zi-0.01z i 2 +(3.186-0.129z i -3.48TO c )TO c ,(l) 
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and which reproduces the data of this work for M c 2; 0.26 
so perfectly that there is no need to establish a new rela- 
tion. Together with Eq. (||) below, Eq. (Q) can be directly 
integrated to obtain L*(r). For the purposes of synthetic 
AGB evolution, the core mass at the instant of time when 
the helium flash occurs can be approximated by 



McRGBtip = (0.454 - 0.023zi + 0.059AL) x 



(2) 



(1 - (0.033 - 0.008zi)e 



MA 



which has a maximum of A/ c .RGBtip approximately at an 
initial mass of 

Mi,max = - ln(0.0169 - 0.005zi) - 2.935 - 0.022^. (3) 

For Mi<M iiinax , Af C:RG Btip = Mc,RGBtip(M iimax ) must be 
used. 

When M CiRG Btip ^ 0.35 according to Eq. (2), he- 
lium ignites under non-degenerate conditions, and there 
is no real RGB. This is the case for M\ <; 2.5. For 
My < 1.5, there is no systematic difference between OPAL 
and LAOL calculations (AAf c < 0.005). The transition 
from degenerate to non-degenerate ignition turns out to 
be relatively model dependent. However, this does not af- 
fect the AGB, since for Mj <; 2 the RGB tip luminos- 
ity does not exceed 1500Lq, and the transition timescales 
from degenerate to non-degenerate are relatively short, so 
that almost no mass is lost. 

The following core helium burning phase does not di- 
rectly influence the AGB, since the luminosities are too 
low to drive significant mass loss, except for blue hor- 
izontal branch stars with an already very small mantle 
mass. These may directly evolve upward in the HRD and 
form AGB manqucc stars, which are not further consid- 
ered here. 

3.2. The E-AGB 

After core helium exhaustion, the helium shell initially 
burns quiescently. The first TP occurs when the He- 
exhausted core approaches the H-shell, at Mh — Mh c 
10~ 2 . L* rises from the relatively low values typical for 
convective core helium burning up to the high AGB lumi- 
nosities. Since the evolutionary time scales on the E-AGB 
are long compared to TP recurrence times, and since for 
stars with AL 4 the maximum luminosity attained on 
the E-AGB exceeds the luminosity at the first TP, a con- 
siderable amount of mass may be lost already on the E- 
AGB. This is interesting with respect to the initial- final 
mass relation, since for M; « 4 there are indications that 
the final mass is less or comparable to the core mass at 
the first full thermal pulse (VW93). 

On the E-AGB, there are competing contributions 
from H- and He-burning, depending on Afj and the evo- 
lutionary phase, so that it makes no sense to use a 
CMLR. Instead, the normalized E-AGB luminosity Ze 



is described directly as a function of the E-AGB phase 
0e = (t - t )/(ti - r ), < 4>e < 1: 

' ^0 



Ie = 



h-l 



ii 



= V(j)E + (1 - V) 



4- 



(4) 



The mass of the helium exhausted core at the beginning 
of the E-AGB, AfHe,o(Ali, Z- 1 ) 1 is calculated from Mne.o = 
max (<7i + (72 x Mi, + 54 x AL), where the coefficients q 
are listed in Table 1, for pop I and II, and can be linearly 
interpolated in \gZ{. In case of "breathing pulses", Mhc,o 
may be larger by up to (0.11 — O.OlMi). 

Once AIh c ,o is known, the duration of the E- 
AGB, T\ — To, and the two parameters v and (3 in 
Eq. (Q) are approximately piecewise linear functions of 
•^He,o alone, independent of Z;. We find that lg(ri — 
r ) = max ((8.4 - 5.0 x Af Hc ,o), (7.6 - 3.1 x M He ,o)), « = 
max (0.19, (-0.23 + 1.35 x M He ,o)) 

and (} = max ((16.3 — 34.1 x Mh c ,o), 5.7), for AfH e ,o vary- 
ing between 0.14 and 0.56. 



Pop 


<?1 


<?2 


<73 


94 


I 


0.168 


0.0174 


-0.175 


0.125 


II 


0.150 


0.0000 


0.050 


0.120 



Table 1. The coefficients to calculate Mh c ,o for Zi — 0.02 and 
10" 4 . 



4. The luminosity on the TP- AGB as a function of 
time 

4-.1. The core mass - luminosity relation (CMLR) 

One of the most important relations that come out of evo- 
lutionary calculations is the maximum luminosity during 
quiescent H-burning, here referred to as La, for which the 
following expression according to Sect. 2.3 is derived: 



L A = (18 160 + 3980^) (M c -0.4468) 

, ^q2.705 + 1.649M c x 

x ^g0.0237(a-1.447)M^ M^(l-o _z 
_ ^q3.529- (M c . -0.4468)AM c /0.01 



(5a) 
(5b) 

(5c) 

(5d) 



The first term is recognized as the usual linear CMLR, 
but the slope is flatter than usual. The term (5b) provides 
a correction: when added to (5a), with (5c) set to unity 
and ignoring (5d), the sum yields for M c ^ 0.6 the same 
numerical values for La as most CMLRs that were fitted 
for low core masses. On the other hand, for Af c 0.95 
this sum approaches the formula given by IT, which is 
applicable for massive cores only (dotted lines in Fig. 0) 



(i 
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The term (5c) corrects for HBB, which is characterized 
by a steep initial increase of La, and a drop when the man- 
tle mass M m is reduced (B95, their figure 9). The calcula- 
tions by DAM96 for Mi < 5 show that HBB operates only 



for massive cores and mantles. The factor (1 — e 



-AM C /0.01 



mimics the initial increase of L*, which levels out after 
about the first ten TPs. HBB strongly depends on the 
MLT parameter, with the results of DAM96 that made 
use of the CMT being described by «mlt = 2.75. The 
data given by VW93 (their figure 12), not used in Eq. (5), 
are remarkably well reproduced by the choice a = 2.25. 

The last term, (5d), corrects for the TOE, which is 
quite important for M c Ss 0.65, and independent of Z. 
Figure | clearly shows that the CMLRs by WZ81 and 
LA86 are influenced by the TOE for low core masses. 

4-2. The core growth equation 

The equation that describes the core growth reads: 



where q, the mass burnt per unit energy release, turns out 
to be slightly metallicity dependent due to the contribu- 
tion from the pp-cycle for low Z. The calculations with a 
nuclear network yield the mean value q — ((1-02 ± 0.04) + 
0.017z)10- n (M /L yr). The luminosity produced by H- 
burning that enters Eq. (||) is smaller than the total lumi- 
nosity, since gravitational energy release due to the core 
shrinking and the flow of burnt mantle material down to 
the core, together with relic helium burning, additionally 
contribute during quiescent H-burning: 



lg(L H ,A/L A ) = -0.012 -10 



-1.25-113AA/ C 



0.0016M m .(7) 



This equation was derived for the standard value of the 
MLT parameter, which implies that the value of La to be 
used in Eq. (7) should be calculated from Eq. (5), with 
the term (5c) put to unity as the standard models did not 
experience any HBB. 

4-3. Further relations for luminosity extrema 

The first observable effect of a TP is that the layers above 
the helium shell source expand, thus the hydrogen shell is 
extinguished and the total luminosity drops. This "rapid 
dip" is the more pronounced the smaller the mantle mass 
is, since the inner layers of the mantle act as a reservoir 
of thermal energy that partly compensates the initial fast 
luminosity drop. The rapid dip is utterly unimportant for 
the secular evolution due to its very short duration, but 
during this phase |dL*/dt| reaches the highest values a 
single low-mass star can achieve. A fraction of roughly 
10~ 3 of all pulsating AGB stars are expected to show cor- 
responding rapid period changes. The fit formula given 
below strongly depends on M m even for low core masses: 
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Fig. 3. The maximum luminosity during quiescent hydrogen 
burning as a function of the core mass and metallicity. Dots: 
data points extracted from the calculations with «mlt = 1-5; 
O: data for «mlt > 1.5 from the calculation with Mi = 5, 
DAM96 and B95. Full lines show the fit to the data given 
by Eq. (5). : inde pende nt data, not used for the fit, from 
Vassiliadis & Wood (1994) for post- AGB nuclei, i. e. objects 



with a vanishing mantle mass. The fit for this case, and without 
TOE, is shown by the dotted lines. The other lines represent 
CMLRs by va rious authors. For the formul ae by Boothroyd 
& Sackmann ( [L988| , BS88) and Lattanzio fll98£|, LA86) the 



dependence on /i, Zcno and Y, respectively, has been taken 
into account approximately. 



L B = 
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" x (8a) 



= (4.81(M C - 0.4865) + M°- 393 / 

/^q2.879 + Af c _ 2Q3.287- (M C-0 -0.4865)AA/ C /0.01 



.(8b) 



After the rapid dip, for M c ^ 0.7 the layers above the in- 
active H-shell start to contract, therefore releasing gravi- 
tational energy; for heavier cores, L* simply follows I/He 
after a short thermal adjustation phase. In either case a 
"rapid peak" emerges, in which the maximum total lumi- 
nosity exceeds the quiescent H-burning luminosity prior to 
the TP, except for the very first weak TPs or in the case of 
strong HBB. Note that in many earlier calculations these 
peaks were suppressed for large mantle masses due to an 
insufficient numerical resolution. It turns out that a HBB 
term does not improve the fit, and is therefore left out: 

L c = 59 200 M c 2 - 10 950 - (9a) 

2Q2.559+ 1.951A/ C , + (44.7+ (2254Z l -147.6)M c ,o)AM c /g^N 

Here Eq. (9b) accounts for the TOE, which vanishes for 
AAf c <; 0.03. In the asymptotic regime, the data nicely 
follow a quadratic relation (Fig. ||). Lq is an interesting 
quantity for two reasons. First, the rapid peak is the only 
mechanism capable of populating a high luminosity tail in 
the luminosity function of a sufficiently large sample with 
Mj < 3. Second, if it should be true that an avalanche 
mass loss starts as soon as some luminosity threshold is 
exceeded (Tuchman et al. 1979), this would take place for 
the first time during a rapid peak, unless HBB is operating 
so effectively that La > Lc- It is a question of the inter- 
play between the mass loss time scale and the duration 
of the peak, if such a gasping mass loss episode influences 
the secular evolution. 

Finally, the helium burning declines and the H-shell 
recovers, giving rise to an extended and long-lasting "slow 
dip". During the quiescent part of a TPC, Ld < I* < La- 
The fit, like Eqs. (8) and (9), needs no HBB correction: 



L D = (76360 + 6460z;)(Af c - 0.3881) 

2q1.55 + 2.11M Ci0 + (59. 7-119. 8M C , )AA/ C 



4-4- The core mass-interpulse time relation 



(10a) 
(10b) 



The second most important quantity to come out of evo- 
lutionary calculations is the core mass-interpulse time re- 
lation. The following relation is found: 

lgrip = (-3.628 + 0.1337zi)(M c - 1.9454) - (11a) 

-^Q-2.080-0.353zi+0.200(A/ m +Q-1.5) _ (lib) 
-^q-0.626 - 70.30(M Oj0 -Zi)AM c (He) 

Again, Eqs. (lib) and (11c) account for HBB and the 
TOE, resp. Interestingly, Eq. (11) predicts the recurrence 
time of pre-pulses with M c < M C: o (crosses in Fig. ||). For 
M c ^ 0.7 and in the asymptotic regime, the results of 
the present work agree well with relations for T lp given by 
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Fig. 4. The luminosity peak after a TP as a function of the 
core mass for i?i = 0.02. Dots: data points from the calculations 
with omlt = 1.5 and used for the fit; O: omlt = 2, : data from 
DAM96, both not used for the fit, which is shown by full lines. 
The appropriate values for M c fi are used for the individual 
sequences. The dotted line depicts La without TOE and HBB 
as in Fig. |§| The dashed line is the asymptotic relation Eq. 
(9a). 
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Fig. 5. The minimum luminosity at the beginning of a TPC 
as a function of M c for Z\ = 0.02. Dashed line: Eq. (10a), other 
symbols: see Fig. m 



BS88 and WZ81. However, the latter disagree for M c > 0.8, 
and give lower and upper limits only, respectively. 

Compared to the asymptotic exponential relation, the 
TOE reduces r; p by almost a factor of two. The influ- 
ence of Z\ on Tjp is rather strong, in the sense that the 
interpulse times increase with decreasing metallicity. This 
also means that the He-shell accretes material processed 
by the H-shell for a longer time, so that the TPs become 
more violent. The reason is that the plasma in a typical 
pop II He-shell is more degenerate than in a pop I-shell. 
An important consequence is that the third dredge-up oc- 
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Fig. 6. The interpulse time as a function of the core mass and 
metallicity. Data points taken from the calculations, used for 
the fit, are shown as dots, the values given by Eq. (11) as full 
lines, o: data from DAM96, and from the run with omlt = 2, 
Mi = 5 (M c <; 0.9). x: recurrence time of pre-pulses, not used 
for the fit. Dashed and dash-dotted lines: relations from BS88 
and WZ81, respectively. Dotted lines: the asymptotic relation 
Eq. (11a). 



curs the earlier, i.e. for lower core masses, the lower the 
metallicity. 

A further quantity needed to construct secular 
lightcurves is the time td that elapses between a TP and 
the slow dip, i.e. when = Ld- We find 

lg(r D /r ip ) = -1.01 + 0.20 * (M c - 1). (12) 

4--5. The core mass at the first thermal pulse 

There is one quantity of outstanding importance left to be 
discussed, the core mass at the first TP, i. e. M Ci o(Mj, Z\). 
It determines the minimum remnant mass, and is there- 
fore closely connected to the initial-final mass relation. For 
M\ 5s 2, the stars develop very similar degenerate He-cores 
(see section 3.1), and all end up with M Cj o ~ 0.55. On 
the other hand, massive stars with Mj <; 4.5 experience 
the second dredge-up, which shifts the H-He-discontinuity 
inward (in mass), with the resulting M Cj o being a linear 
function of M\. In between, there is a transition region, 
which is rather strongly influenced by computational de- 
tails. For stars in this transition region, dM Cj o/dMi <^ 0.25, 
and since increases strongly with M c , and the mass loss 
even much more steeply with L*, such transition stars, 
despite their larger mass reservoir, have a shorter AGB 
lifetime than slightly less massive counterparts. In WW94 
the following parameterization was introduced: 



M C) o(M0 

/(MO = (l + o 



(-p 1 (M i -p 2 ) 2 +p 3 )/ + 
{p 4 Mi + p 5 ) (1 - /), where 

Mj-pgN 



(13a) 
(13b) 

(13c) 



Equation (13a) mimics M Cj o for stars with Mi < 2.5 that 
experience the core helium flash, which is almost constant. 
Eq. (13b) accounts for the linear ingress of the second 
dredge-up for Mj ^ 4.5, and Eq. (13c) mediates. All pa- 
rameters can be interpolated linearly in lg-Zi, hence in 
Table |^ their numerical values are given for two standard 
metallicities only. 

A compilation of data from various authors is shown 
in Fig. 0. There is a tendency that more recent opacities 
yield higher values of M Cj o- Blocker's (B95) data points 
lie systematically below all other's. Obviously there is no 
systematic effect due to the inclusion of semiconvection 
(LA86, DAM96). The compilation by GdJ systematically 
deviates from the data for the LMC metallicity in the 
sense that for low-mass stars it predicts too large values 
of M c n. 
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0.0294 


1.478 


0.550 


0.0634 


0.572 


3.193 


0.260 


II 


0.0213 


2.589 


0.592 


0.0324 


0.790 


2.867 


0.260 



Table 2. The coefficients of Eq. (13) for Zi = 0.02 and 10" 
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4-6. Luminosity variations during the thermal pulse cycle 
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Fig. 7. The core mass at the first TP as a function of the 
initial mass and metallicity. The symbols show data from vari- 
ous authors (CC90: Castellani et al. 1990). The lines show the 
relation 13a (dots) and the one used by GdJ (dashed). 



In previous synthetic AGB calculations, the luminosity 
variations during a thermal pulse cycle were either ne- 
glected (Rcnzini & Voli 1981), or described by simple step- 
functions (GdJ, Marigo et al. 1996). This is clearly insuffi- 
cient if accurate theoretical luminosity functions are to be 
predicted. This section completes the description of -L*(r) 
on the AGB. 

At the instant of time when L„ = L^>, i. e. in the slow 
dip, the H-shell is recovering from the preceding TP, and 
in ~ £hc- It turns out that the luminosity produced by 
H-burning that enters Eq. (^J), when normalized to the 
value predicted by the CMLR Eq. (5) and Eq. (0), has an 
almost universal form (see also WZ81), independent of Z{. 
Let 4> ee r/T ip (4> = at the TP initiating the TPC, when 
M c = M c tp). Then we define 



M0) = 



i 



-M 2 



J H,CMLR 



(M c (0),M C)O ,M* 



1 



-0H 



(14) 



where /?h = 93.11(M c .tp — 0.4569). There are small sys- 
tematic deviations for the first TP of each sequence, which, 
however, have no influence on the integration of Eq. (|6|). 

depends similarly on cj). In the following, the ar- 
guments of quantities given by a CMLR are omitted; the 
index "TP" denotes a quantity as given by the correspond- 
ing CMLR for M c = M c tp- The form function h^ is given 

by 



M0) 



M0) 



^A,TP 
L A (Mc(0)) 



^D,TP 



1 



-mi 



La 



TP 



^D.TP 



=.-/3l 



(15) 



where /3 L = 7.95 and b L = 2.13 (1.331 - M c , T p), and ee 
max{(/)— 0.1, 0} accounts for the behaviour shortly after the 
slow dip. Note that h L (l) > 0.999« 1, and that 0.1<</><1 
due to Eq. (|l|). 

Let now ee (t — tc)/t iv {M Ci M Cj0 , M„) be a shifted 
TPC phase, so that <j> — at the rapid peak when = 
Lq. The peak itself can be approximated by a parabola, 
followed by an exponential decline. The luminosity from 
the rapid peak until L* = Ld can well be described by: 



'aft, 



(16) 



hi 

4a ' 



if </>< 



A N (l-/ N )-/ N (-60+£] else. 



4a 



An 
b 



The physical meaning of An and /n is, that when the lu- 
minosity Z» = lg Lq — An is reached, the time scale for 
the decline changes from T lp /b to Ti p /(f^b). The parame- 
ters lga,lg6, lg An and lg/N are functions of Zi, M c and 
M Ct Q and can again be described by linear relations with a 
TOE correction for lg a and lg b. They obey the simplified 
relations 



lga = 7.23- 4.40 M c , 



(17) 
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Igb = 3.30- 2.91 M c , 
lg An = 0.08- 1.50 M c , and 
lg/ N = -1.84+ 1.79 M c . 



(18) 
(19) 
(20) 



Note that the duration of the rapid peak does not scale 
like Tip. Relative to Ti p , it lasts the longer the more massive 
the core is, i. e. the shorter ri p . 

5. Properties of the mantle 

5.1. The method 

Several physical aspects of AGB evolution depend on 
properties that are almost exclusively determined by the 
mantle. The justification is that as long as there is a ra- 
diative transition region between the core and the lower 
boundary of the mantle convection zone, the influence of 
the two outer boundary conditions on the core is negligible 
(WZ81). On the other hand, since the luminosity is essen- 
tially constant, the core, which prescribes L, can to first 
order be considered to be just a gravitating point source 
for the mantle. 

Now consider L*,Af*,M c and the chemical composi- 
tion to be given. The latter, owing to mixing during earlier 
evolutionary phases, is homogeneous up to and including 
the upper parts of the H-shell. Two outer boundary con- 
ditions for the mantle are given immediately: the effective 
stellar surface radius P* is given by the Stefan-Boltzmann 
law and from some model for the optically thin outer en- 
velope one can derive P*(P*, M», T„). By assuming a value 
for T*, the four stellar structure equations for L,T,P (or 
p) and R as functions of the Lagrangean mass coordi- 
nate M r can now be integrated inward. The transition to 
the core is well defined (see also Wagenhuber & Tuchman 
1996| ) by a very steep rise of T and P by one and four, 



resp., orders of magnitude or more. In this way, a relation 
M C (T*) is defined that can be inverted numerically to yield 
T* as a function of all other quantities, defining a three 
dimensional manifold in the four dimensional parameter 
space (T„, L„, M*, M c ). Full stellar evolution calculations 
of course establish a connection between L* and M c . In 
order to study pure mantle properties, however, it is much 
more convenient to treat them as independent variables to 
get rid of the complicated time dependent problem one is 
confronted with in full calculations. 

In the following, a set of data relevant for HBB is dis- 
cussed, together with consistent effective temperatures, for 
various model assumptions that made use of OPAL (see 
Sect. |2.2| ) and more recently published opacity tables. 

5.2. The effective temperature 

Figure H illustrates how T*(i*,M*,M c ) depends on the 
various variables and physical model assumptions. Let us 
first discuss the standard case (pop I, M* = 1, MLT; up- 
per left box). For a given core mass, decreases almost 



linearly for increasing I* as known from the Hayashi the- 
ory, the exact value of the slope depending on the opacity 
as a function of P and T, here dt*/dl* ~ —0.2. When Z* 
approaches the Eddington luminosity (pEdd)j t* becomes 
locally almost independent of Z* , and there is even a bend- 
ing towards hotter t* for large M c . Such mantles are no 
longer fully convective. It is obvious that the pressure at 
the photosphere has a minor influence (dotted lines, up- 
per left box), since on the AGB for I* > 3, < 3 for all 
atmosphere models considered. 

If the mantle mass is reduced, i. e. for growing M c with 
M* kept constant, increases. The influence of M c on t* 
is the less pronounced the larger M* is; for M* ^ 2 it 
is already almost negligible. Increasing M*, while M c = 
const, leads to higher effective temperatures and vice versa 
without changing dt*/dk significantly (upper right box). 
This effect is responsible for the trend towards lower effec- 
tive temperatures on the upper AGB as the objects lose 
mass, as long as M m <; 0.2. If M m becomes smaller, the 
former effect takes over, and for M c ss const a shrinking 
total mass leads to higher T*, and thus to the departure 
of proto-PN objects from the AGB. 

Increasing the MLT parameter has nearly the same 
effect as reducing the metallicity (middle row): it shifts 
the whole set of curves towards higher effective tempera- 
tures without changing the influence of M c . It turns out 
that there is no qualitative difference between the MLT 



and the Canuto-Mazzitclli theory (1991, CMT, lower row). 
Quantitatively, there is no single value of qcmt which re- 
produces £*,mlt for all M c and U to within 300K. This 
difference, however, is much smaller than the errors of the 
available observational determinations of t*. Recent data 



(van Leeuwen et al., 1997 ) that make use of hipparcos 
parallaxes, shown in the lower left box of Fig. ^, support 
the lower IVscale. E. g. \ Cyg can well be explained to 
be a low-mass object (M* « 0.75) with a low-mass pro- 
genitor, so that M C: o ~ 0.535, in one of its first TPCs 
(M c = 0.540). 

5.3. The bottom of the convective mantle 

With respect to HBB, T b (L*,M*,M c ) together with p b 
determines the rate of nuclear reactions (see e.g. Scalo et 
al. 1975| ). A fact that has been widely disregarded in the 
past is that also the steepness of the T- and P-profilcs 
is decisive for HBB, since it, together with the mixing 
time scale, determines the amount of matter subject to 
burning. There are two extreme possibilities: On the one 
hand, the fuel could be burnt locally at the bottom of the 
mantle convection zone. On the other hand, if the mixing 
is sufficiently effective, the mantle material as a whole at 
any instant of time could be subject to nuclear processing. 
The local mixing time scales, which are very poorly known 
from the MLT, determine the actual situation between 
these extremes. As long as the uncertainties due to the 
convection formalism dominate to such an extent, it is 
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Fig. 8. The effective temperature as a function of L*, M*, M c and Z, and the influence of various model assumptions by means 
of the example M„ = 1. The standard case is depicted in the upper left box; the full lines are T*(L*) for M c — 0.5, . . . , O.95M . 
This case is repeated in all other boxes as dotted lines for comparison. In the upper left box itself, the dotted lines are obtained 
by reducing the pressure at the outer boundary condition by a factor of three as compared to the usual relation P*(L*, M*,T»). 
In the lower left box, recent observational data (see text) are overplotted with la error boxes. 
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Fig. 9. The temperature and some quantities characterizing the overall conditions at the bottom of the convective mantle 
(see text) as functions of L*,M* and M c for the standard case (pop I, MLT, a — 1.5, left column) and for various physical 
assumptions (right column). Left: the thick lines show the data for M c = 0.5 and different total masses as indicated. The 
thin lines illustrate the effect of M c (for M* = 1, the highest core mass is 0.9Mq, in all other cases I.IMq). Right: The same 
quantities for M* = 7.5 and M c = 0.7 and 1.1 (left resp. upper set of curves) for various physical model assumptions. 
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well justified to use parametrized descriptions instead of 
time consuming full stellar evolution calculations. 

The most important quantity is of course the temper- 
ature at the base of the convective mantle (Fig. |9|, upper 
left box). The present work qualitatively confirms the re- 
sults obtained by Boothroyd et al. (1993) by means of 
full stellar evolution calculations, and work out the ba- 
sic dependencies more clearly, like the shift of Tb caused 
by breathing pulses (their Figure lc), which can be easily 
explained as the effect of an increased core mass. 

For sufficiently large total masses M» <; 2, % is roughly 
independent of M*. Cool bases (<b 7) occur for I* 4. 
If the luminosity increases beyond some critical value de- 
pending on M c , Xb either rises approximately like Tb oc L\ 
for low core masses, or drops, since the inner parts of the 
mantle become radiative for large M c and small M* (in 
Fig. |, this is the case for U > 4.2, M* = 3 and M c > 0.7). 
At tb ~ 7.9, the rise flattens, unless U gets very close to 
the Eddington luminosity. For a given set of opacity tables, 
"hot bottoms" can occur for sufficiently large luminosities 
only if the ratio |p is less than some constant that de- 
pends solely on Z and the convection description. A useful 
approximation is 



Mo < f 0.2 Q!MLT ' 

M* \0.5a C MT 



5Z, for MLT, and 

- hZ - 0.08, for CMT. 



(21) 



The influence of the physical models is only moderate (Fig. 
^, upper right box) and leads predominantly to a shift of 
tb(J*) along the Z*-axis, which is equivalent to a rescaling of 
M c . Every change of the physical model input that leads 
to higher effective temperatures, like lower opacities at 
low temperatures, lowers the threshold for Z* above which 
HBB takes place. Changes of the opacities k(T, P) can be 
compensated by varying a appropriately, as can be clearly 
seen for the pair of cases shown with long dashed and 
thick dotted lines. For the former, CMT and «cmt = 0-80 
was used, for the latter MLT, a = 1.5 and recent OPAL 
opacities with low temperature extensions according to 



Alexander & Ferguson (1994, labelled AF94 in the figure). 
For the pop II models, which imply a drastic reduction 
of k at low temperatures, also the plateu value of tb is 
increased significantly (dotted lines in Fig. [| upper right 
box). 

In order to calculate nuclear reaction rates, the den- 
sity at the convective base must be given, too. Since the 
mantle convection zones are almost entirely adiabatic, pb 
also fixes the entropy. For a sufficiently large total mass, 
and unless L* approaches the Eddington limit, the simple 
relation 



T 3 
lg- b 



PhL* 



18.9 - 0.4M C 



(22) 



holds remarkably well over one order of magnitude of L* 
(Fig. ||, upper middle row) and may be used to estimate 
pb- 



As has been already mentioned, the quantity 4rfr~ 

(Fig. ^ lower middle row) determines the amount of mat 
ter subject to nuclear processing. It increases approxi- 
mately as oc L 3 , and even superexponentially when L* — > 
iEdd- This means, that the effect of reaching very high 
base temperatures could be entirely compensated by the 
fact that almost no fuel is left. It must be emphasized, 
that the nuclear yields will depend extremely sensitively 
on the mixing time scales in such a situation, which are 
not correctly predicted by the MLT! 

In the lower row of the figure, the adiabatic gradient 



V 



ad.b 



= SL„ is shown. Let h= ^ 



with 1.8 <; 



dlnPlad,b — d In T 

h < 3, denote the steepness of the density profile. Then 
V a d.b, the ratio j3 of the gas pressure to the total pressure, 
including that of radiation, and h are (exactly) connected 
by the relations 



2h-7 



ad,b 



3/3 



4 2h 2 -6/i-3 16 - 12/3 - 



I/? 2 



(23) 



as long as degeneracy plays no role for the EOS. In any 
case, pb < 10g cm~ 3 , so that this assumption is well 
fulfilled. The uniform decrease of V a d,b as increases 
demonstrates the increasing contribution of the radiation 
pressure to Pb- For L* — > i-Edd, /3 — > and V a d,b - * 0.25. 
ib('*) reaches the plateau value of % w 7.9 approxi- 
mately when the radiation pressure dominates, /3 < 1/3, 
or V ad ,b < 0.255. 

6. Summary 

We present analytical formulae that describe the evolution 
of a star on the AGB. This outdates all previous work on 
this subject up to now for several reasons. First, all formu- 
lae are based on a (nearly) homogeneous set of up-to-date 
full stellar evolution models that cover all relevant masses 
and a large range in metallicity. Second, the analytical for- 
mulae take into account several important features found 
in the full calculations, and largely neglected so far in pre- 
vious analytical descriptions: the secular variations of the 
luminosity during a thermal pulse cycle, turn-on effects 
during the first few pulses, and hot bottom burning. What 
these descriptions and formulae can possibly not account 
for, is the influence or feed-back of dredge-up (or nucle- 
osynthesis in general) on the evolution of the star. 

In future work we plan to include a description of the vari- 
ations of the various chemical species during 1st, 2nd and 
3rd dredge- up and HBB, which were not addressed in the 
full stellar evolution models. 

We are currently writing a numerical code to implement 
the analytical formulae described in this paper. When the 
chemical description is included we are in a position to 
improve upon the AGB population synthesis models of 
GdJ or Marigo et al. (1996). With such a code we are 
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also in a position to investigate new data on AGB stars in 
extragalactic systems that will become available with the 
current and new generation of 8-10m telescopes. 
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